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LONG-TERM  GOALS 

The  ultimate  goal  is  to  develop  direct  simulation/physics-based  forward  and  inverse  capabilities  for 
radiance  prediction  in  a  dynamic  ocean  environment.  This  direct  simulation-based  model  will  include 
and  integrate  all  of  the  relevant  dynamical  processes  in  the  upper  ocean  surface  boundary  layer  into  a 
physics-based  computational  prediction  capability  for  the  time-dependent  radiative  transport. 

OBJECTIVES 

To  include  and  integrate  relevant  dynamical  processes  in  the  upper  ocean  surface  boundary  layer 
(SBL)  into  a  physics-based  computational  prediction  and  inverse  capability  for  the  time-dependent 
radiative  transport: 

•  Develop  direct  simulation  of  upper  ocean  hydrodynamic  processes  and  forward  prediction  of 
radiative  transfer 

•  Obtain  understanding,  modeling  and  parameterizations  of  dependencies  of  oceanic  radiance  on 
the  surface  wave  environment 

•  Provide  guidance  for  field  measurements  and  obtain  cross  validations  and  calibrations  with 
direct  simulations  and  modeling 

•  Provide  a  framework  for  inverse  modeling  and  reconstruction  of  ocean  surface  and  above  water 
features  based  on  sensed  underwater  radiance  data 

To  reach  these  objectives,  we  had  and  would  continue  to  have  a  close  collaboration  with  Professor 
Lian  Shen  of  the  Johns  Hopkins  University  (JHU)  on  the  modeling  of  free  surface  turbulence 
roughness. 

APPROACH 
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A  simulation  approach,  based  on  direct  physics-based  simulations  and  modeling,  is  applied  to  solve  the 
problem  of  ocean  radiance  transport  (RT)  in  a  dynamic  ocean  SBL  environment  that  includes 
nonlinear  capillary-gravity  waves  (CGW),  free-surface  turbulence  (FST)  roughness,  wave  breaking, 
and  bubble  generation  and  transport.  The  radiative  transport  is  governed  by  Snell’s  law  and  Fresnel 
transmission  at  the  water  surface,  and  absorption  and  multiple  scattering  in  the  water  underneath.  The 
complex  dynamic  processes  of  the  ocean  SBL,  the  nonlinear  CGW  interactions,  the  development  and 
transport  of  FST,  and  the  generation  and  transport  of  bubbles  are  modeled  using  physics-based 
computations.  The  modeling  of  these  hydrodynamic  processes  is  coupled  with  the  computation  of 
radiative  transport. 

Nonlinear  CGW:  An  efficient  phase-resolved  computational  approach  based  on  Euler  equations  is 
used  to  compute  spatial  and  temporal  nonlinear  evolution  of  capillary  and  gravity  waves.  This 
computational  tool  builds  on  an  efficient  high-order  spectral  method  that  we  developed  for  direct 
simulations  of  nonlinear  gravity  wavefield  evolution.  Nonlinear  gravity-gravity  and  gravity¬ 
capillary  wave  interactions  are  accounted  for  up  to  an  arbitrary  order  in  the  wave  steepness.  With 
high  computational  efficiency,  this  approach  enables  phase-resolved  simulations  of  large-scale 
nonlinear  CGW. 

FST-Wave  Interactions :  Navier-Stokes  equations  based  DNS  is  employed  to  resolve  all  eddies  in 
free  surface  turbulence.  LES  and  LWS  are  used  to  compute  large  eddy  and  large  wave  components 
explicitly,  with  effects  from  small-scale  motions  being  represented  by  subgrid-scale  (SGS)  models. 
Fully  nonlinear  viscous  free-surface  boundary  conditions  are  imposed.  Effects  of  surfactants  are 
captured  through  the  Plateau-Marangoni-Gibbs  effect  with  surfactant  transport  directly  simulated. 

Steep  and  Breaking  Waves:  A  Navier-Stokes  equations  solver  for  fully  coupled  air-wave 
interactions  with  a  level-set  method  for  free  surface  tracking  is  employed  to  compute  the  details  of 
free  surface  signature  and  dissipation  due  to  steep  and  breaking  waves. 

Bubble  Transport  in  CGJV  and  FST:  Direct  simulation  is  developed  to  compute  bubble  motion  in 
CGW  and  FST  environment.  Both  Lagrangian  and  Eulerian  approaches  are  used  to  trace  bubble 
trajectories.  Bubble  motion  is  subject  to  forces  due  to  added  mass,  buoyancy,  drag,  lift,  and  fluid 
stress  gradients  arising  from  the  continuous-phase  acceleration.  Bubble  source  is  determined  based 
on  experimental  measurements  and/or  existing  data. 

Radiative  Transfer^  in  CGW  and  FST:  Monte  Carlo  simulation  of  radiance  transfer  (RT)  is 
developed  with  the  free  surface  deformation  obtained  from  direct  CGW  and  FST  computations. 
The  effects  of  absorption  and  multiple  scattering  on  RT  are  included.  Bubble  scattering  effects  are 
also  considered  based  on  bubble  distribution  and  transport  and  Mie  theory. 

WORK  COMPLETED 

•  Development  of  nonlinear  CGW  simulation  capability :  Developed  a  direct  simulation 
capability  for  nonlinear  CGW  evolution  by  extending  the  high-order  spectral  method  for 
gravity  waves  to  general  broadband  nonlinear  wavefield  interactions  involving  swell,  seas,  and 
capillary  waves.  The  focus  was  on  the  inclusion  of  capillary  waves,  long-short  wave 
interactions,  wave  breaking  dissipation,  dissipation  due  to  free  surface  turbulence  and  energy 
input  by  wind. 
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•  Development  and  validation  of  unpolarized/polarized  Monte  Carlo  RT  simulation : 
Developed  a  three-dimensional  coupled  atmosphere-ocean  Monte  Carlo  radiative  transfer 
simulation  capability  for  both  polarized  and  unpolarized  lights.  The  RT  simulation  is  time 
independent,  but  accounts  for  the  effect  of  complex  unsteady  three-dimensional  ocean  surface 
including  CGW  and  FST  roughness.  In  RT,  multiple  refractions  at  ocean  surface,  total  internal 
reflection,  all  orders  of  multiple  scattering,  and  scattering  and  absorption  of  both  water 
molecules  and  marine  aerosols  are  all  considered.  In  order  for  practical  applications,  various 
techniques  including  the  use  of  biased  sampling  algorithms  and  parallelization  of  the  code 
(with  MPI)  were  employed  to  speed  up  the  program.  The  Mont  Carlo  RT  simulations  were 
validated  by  comparisons  to  available  experimental  data  and  other  model  predictions. 

•  Investigation  of  radiance/irradiance  and  polarization  in  CGW/FST  environment  and  under 
linear  and  nonlinear  waves.  Investigated  the  patterns  of  underwater  radiance  and  polarization 
induced  by  the  CGW/FST  dynamic  ocean  surfaces  and  their  relations  to  solar  conditions  and 
ocean  IOPs.  Both  linear  and  nonlinear  surface  waves  were  compared  in  the  RT  simulations  in 
order  to  understand  the  effects  of  nonlinearity  of  ocean  surface  on  the  underwater  light  fields. 
For  irregular  linear  and  nonlinear  ocean  surface  waves,  the  spatial  variation  of  downwelling 
irradiance  was  studied. 

•  Development  of  analytical  GP  model  of  probability  density  function  (PDF)  of  underwater 
downwelling  irradiance.  In  order  to  understand  variability  of  the  fluctuating  underwater 
downwelling  irradiance  due  to  dynamic  ocean  surface  focusing  effects,  we  focused  on  the 
probability  density  function  (PDF)  of  the  downwelling  irradiance  and  developed  an  analytical 
Gauss-Poisson  (GP)  model  to  quantify  the  PDF  at  the  arbitrary  depths.  The  model  predictions 
were  compared  to  experimental  data  and  showed  a  good  agreement  with  data. 

•  Development  of  3D  Monte  Carlo  code  of  radiative  transfer  for  variable  refractive  index  and 
primary  application  to  ocean  turbulence.  Studied  the  unpolarized  light  radiative  transfer 
properties  within  an  ocean  turbulent  flow  by  applying  the  inhomogeneous  radiative  transfer 
equation  with  variable  refractive  index.  A  hybrid  3D  numerical  program  combining  Monte 
Carlo  method  and  ray-tracing  technique  was  developed  to  simulate  the  RT  process  in  the 
inhomogeneous  ocean  turbulent  flow. 

RESULTS 

Validations  of  the  MC  RT  simulation  capability:  We  developed  a  3D  Monte  Carlo  (MC)  RT 
simulation  capability  for  unpolarized  and  polarized  light  for  the  atmosphere-ocean  system.  The  model 
has  been  systematically  validated  by  direct  comparisons  with  existing  theories  and  numerical  model 
predictions  and  with  the  RaDyo  field  measurements.  Figure  la  shows  the  comparison  of  unpolarized 
radiance  distribution  between  MC  RT  simulation  and  the  prediction  by  the  invariant  imbedding 
method;  and  figure  lb  shows  the  comparison  of  degree  of  polarization  between  MC  RT  simulation  and 
RaDyO  Hawaii  experimental  data  at  the  depth  z=-8  m.  Both  model-to-model  comparison  and  model- 
to-data  comparison  show  good  agreements  between  MC  predictions  with  other  models  and 
experimental  measurements. 

Effects  of  ocean  surface  waves  on  underwater  radiance  and  polarization  pattern:  Using  the  MC 
simulation  capability,  effects  of  the  ocean  surface  roughness  on  the  underwater  light  polarized  light 
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fields  were  intensively  examined.  Figure  2a  shows  the  dependence  of  the  ratio  of  maximum  mean 
value  of  degree  of  polarization  to  the  sky  degree  of  polarization  <P>max/Psky  on  the  mean  square  slope 
(MSS)  of  the  ocean  surface  waves.  Four  different  depths  are  considered  and  they  are  z  =  -0.3m,  -2.2m, 
-4.9m,  and  -7.5m.  The  solar  incidence  is  0SUn=65o  and  (pSun=0o;  the  light  wavelength  k=532  nm;  IOPs 
were  obtained  from  the  Santa  Barbara  Channel  (SBC)  experiment  which  are  total  beam  attenuation 
coefficient  c=0.610  m'1  and  total  single  scattering  albedo  eo=0.873;  Petzold  phase  function  was  used. 
It  can  be  seen  that  increasing  surface  roughness  leads  to  smaller  degree  of  polarization  and  the  relation 
between  the  maximum  underwater  degree  of  polarization  and  MSS  are  approximate  linear.  The 
dependence  of  maximum  degree  of  polarization  becomes  weaker  at  deeper  locations,  indicating  that 
the  effect  of  ocean  surface  waves  on  polarization  diminishes  with  detector  depth. 

Consideration  of  nonlinearity  of  surface  waves  is  of  critical  importance  in  the  study  of  underwater 
irradiance  as  it  affects  the  profile  of  ocean-atmosphere  interface.  Figure  2b  compares  the  effect  of 
linear  and  nonlinear  waves  with  same  MSS  and  different  MSS  respectively  on  polarization  in  variance 
of  the  mean  value  of  degree  polarization  <P>  at  the  local  maximum  inside  the  Snell’s  window,  at  the 
position  of  0=40°  and  cp=l  10°.  It  is  seen  that  the  ratio  <P>max/PSky  for  the  linear  and  nonlinear  waves 
with  the  same  MSS(~0.04)  are  equal  to  each  other  at  arbitrary  depth.  <P>max/PSky  induced  by  the  linear 
waves  with  MSS=0.03  is  a  little  larger  than  that  of  nonlinear  waves.  The  variance  of  P/<P>  at  the 
local  maximum  point  with  nonlinear  waves  is  larger  than  the  variances  of  linear  waves  with  the  same 
and  smaller  MSS.  Therefore,  it  can  be  expected  that  in  the  case  of  rougher  sea  surfaces  where 
nonlinearity  becomes  more  dominant,  the  variance  of  the  in-water  polarization  is  larger  than  the 
previous  calculation  based  on  a  linear  wave  assumption. 

Variability  of  underwater  irradiance  induced  by  ocean  surface  wave:  To  understand  the  variability 
and  fluctuation  of  the  in-water  light  fields,  we  developed  a  statistical  GP  model  that  analytically 
quantifies  the  PDF  of  the  downwelling  irradiance  under  random  ocean  waves.  The  model  assumes  that 
ocean  surface  as  independent  and  identically  distributed  flat  facets.  The  theoretical  model  captures  the 
characteristics  of  the  PDF  from  skewed  to  near-Gaussian  shape  as  the  depth  increases  from  shallow  to 
deep  water.  It  obtains  a  closed-form  asymptotic  for  the  probability  that  diminishes  at  a  rate  between 
exponential  and  Gaussian  with  increasing  extreme  values.  Figure  3a  shows  the  comparison  the  PDF  of 
normalized  the  downwelling  irradiance  Ed/<Ed>  at  the  depth  z=-2.85m  with  experimental  data  and 
predictions  by  MC  simulation.  It  shows  that  a  very  good  agreement  with  the  three  methodologies  has 
been  reached.  Figure  3b  demonstrates  one  application  of  GP  model  in  understanding  the  how  wind- 
driven  ocean  surface  wave  spectrum  influences  the  PDF  of  underwater  irradiance.  It  can  been  seen 
that  increasing  wind  speed  results  in  the  monotonic  decrease  of  the  normalized  irradiance  standard 
deviation  ax  and  as  increasing  depth  cx  rapidly  approaches  the  deep  water  limit  as  the  downwelling 
irradiance  is  dominated  by  volume  scattering  effect. 

Effects  of  ocean  turbulence  flow  on  underwater  radiance/irradiance  distributions:  We  considered 
the  radiative  transfer  process  within  the  inhomogeneous  ocean  turbulent  flow.  To  solve  the  non- 
homogeneous  radiative  transfer  equation  with  the  variable  refractive  index,  we  developed  the  hybrid 
Monte  Carlo  method  and  ray-tracying  technique  to  account  for  the  light  beam  bending.  Figure  4a 
shows  the  downwelling  irradiance  pattern  at  the  depth  of  z=-80m  as  the  ocean  surface  is  flat.  The 
irradiance  pattern  reflects  the  structure  of  the  ocean  turbulence.  Figure  4b  shows  the  downwelling 
irradiance  pattern  at  the  depth  of  z=-80m  caused  by  the  progressive  ocean  surface  waves  and  the 
turbulent  flow  beneath  the  surface.  It  can  be  seen  that  basic  pattern  of  the  downwelling  irradiance  is 
determined  by  focusing  of  light  at  the  ocean  surface  and  the  pattern  is  blurred  by  the  turbulent  flow. 
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IMPACT/APPLICATIONS 


The  capability  of  accurate  prediction  of  the  irradiance  transfer  across  ocean  surface  and  in  the  water 
may  enable  the  development  of  a  novel  approach  for  accurate  measurements  of  complex  ocean 
boundary  layer  processes  and  reliable  detection  of  the  presence  of  structures/objects  on  or  above  ocean 
surface  based  on  sensed  underwater  irradiance  data. 
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Figure  1:  (a)  Comparison  of  unpolarized  radiance  distribution  between  present  MC  RT  simulation  and  the 
prediction  by  the  invariant  imbedding  method;  and  (b)  comparison  of  polarization  degree  between  MC  RT 
simulation  (black  solid  line)  and  RaDyO  Hawaii  experimental  data  (red  circle). 
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Figure  2.  (a)  Dependence  of  ratio  of  maximum  mean  value  of  degree  polarization  to  sky  degree  of 
polarization  <P>max/Psky  within  the  SnelVs  window  on  the  MSS  of  the  ocean  surface  at  different  depths ,  and 
(b)effects  of  the  nonlinearity  of  ocean  waves:  the  variance  of  the  normalized  degree  of  polarization  at 
different  depths  under  linear  and  nonlinear  waves  with  same  MSS  and  different  MSS. 
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Figure  3  .(a)  Comparison  of  the  GP  theoretical  model  (—)  for  the  PDF  of  the  normalized  downwelling 

irradiance  /=E/<E>  with  Monte  Carlo  (MC)  simulation  (•••)  and  experimental  data  ( - )  at  depth  of 

2.85m ,  and  (b)  wind  speed  effect  on  standard  deviation  of  normalized  downelling  irradiance  Gxfor  different 
depths:  Dcw0=  0.524(—);  1.04  (( - -);  1. 74((- •  -•)  and DCw0  — 


Figure  4.  (a)  Instantaneous  downwelling  irradiance  Ed  at  z=-50  m  within  a  turbulent  flow  under  the  flat 
ocean  surface ,  and  (b)  instantaneous  downwelling  irradiance  Ed  at  z=-50  m  within  a  turbulent  flow  under 
ocean  surface  waves  of peak  frequency  steepness  kpap=0.1  propagating  in  the  same  direction  as  the  shear 

flow. 
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